Deriving large-scale glacier velocities from a complete satellite archive : Application to the Pamir-Karakoram-Himalaya

International audienceMountain glaciers are pertinent indicators of climate change and their dynamics, in particular surface velocity change, is an essential climate variable. In order to retrieve the climatic signature from surface velocity, large-scale study of temporal trends spanning multiple decades is required. Satellite image feature-tracking has been successfully used to derive mountain glacier surface velocities, but most studies rely on manually selected pairs of images, which is not adequate for large datasets. In this paper, we propose a processing strategy to exploit complete satellite archives in a semi-automated way in order to derive robust and spatially complete glacier velocities and their uncertainties on a large spatial scale. In this approach, all available pairs within a defined time span are analysed, preprocessed to improve image quality and features are tracked to produce a velocity stack; the final velocity is obtained by selecting measures from the stack with the statistically higher level of confidence. This approach allows to compute statistical uncertainty level associated with each measured image pixel.This strategy is applied to 1536 pairs of Landsat 5 and 7 images covering the 3000 km long Pamir–Karakoram–Himalaya range for the period of 1999–2001 to produce glacier annual velocity fields. We obtain a velocity estimate for 76,000 km2 or 92% of the glacierized areas of this region. We then discuss the impact of coregistration errors and variability of glacier flow on the final velocity. The median 95% confidence interval ranges from 2.0 m/year on the average in stable areas and 4.4 m/year on the average over glaciers with variability related to data density, surface conditions and strain rate. These performances highlight the benefits of processing of a complete satellite archive to produce glacier velocity fields and to analyse glacier dynamics at regional scales

Ice loss in High Mountain Asia and the Gulf of Alaska observed by CryoSat-2 swath altimetry between 2010 and 2019

We report the first application of a novel approach to retrieve spatially-resolved elevation change from radar altimetry over entire mountain glaciers areas. We apply interferometric swath altimetry to CryoSat-2 data acquired between July 2010 and July 2019 over High Mountain Asia (HMA) and in the Gulf of Alaska (GoA). We bin swath elevation data into 100 x 100 km bins, remove the topography with a reference DEM and generate linear rates of elevation changes for each bin individually using a weighted regression model. We exclude solutions that that did not fulfil a set of quality criteria based on elevation change uncertainties, temporal completeness, interannual changes and stability of regression results. To extrapolate missing data, hypsometric averaging is applied. We find that during the study period, HMA and GoA have lost an average of –28.0 ± 3.0 Gt yr–1 (–0.29 ± 0.03 m w.e. yr–1) and –76.3 ± 5.7 Gt yr–1 (–0.89 ± 0.07 m w.e. yr–1) respectively. Glacier thinning is ubiquitous except for the Karakoram-Kunlun region experiencing stable or slightly positive mass balanc

Spatially and temporally resolved ice loss in High Mountain Asia and the Gulf of Alaska observed by CryoSat-2 swath altimetry between 2010 and 2019

International audienceGlaciers are currently the largest contributor to sea level rise after ocean thermal expansion, contributing ∼ 30 % to the sea level budget. Global monitoring of these regions remains a challenging task since global estimates rely on a variety of observations and models to achieve the required spatial and temporal coverage, and significant differences remain between current estimates. Here we report the first application of a novel approach to retrieve spatially resolved elevation and mass change from radar altimetry over entire mountain glaciers areas. We apply interferometric swath altimetry to CryoSat-2 data acquired between 2010 and 2019 over High Mountain Asia (HMA) and in the Gulf of Alaska (GoA). In addition, we exploit CryoSat's monthly temporal repeat to reveal seasonal and multiannual variation in rates of glaciers' thinning at unprecedented spatial detail. We find that during this period, HMA and GoA have lost an average of −28.0 ± 3.0 Gt yr−1 (−0.29 ± 0.03 m w.e. yr−1) and −76.3 ± 5.7 Gt yr−1 (−0.89 ± 0.07 m w.e. yr−1), respectively, corresponding to a contribution to sea level rise of 0.078 ± 0.008 mm yr−1 (0.051 ± 0.006 mm yr−1 from exorheic basins) and 0.211 ± 0.016 mm yr−1. The cumulative loss during the 9-year period is equivalent to 4.2 % and 4.3 % of the ice volume, respectively, for HMA and GoA. Glacier thinning is ubiquitous except for in the Karakoram-Kunlun region, which experiences stable or slightly positive mass balance. In the GoA region, the intensity of thinning varies spatially and temporally, with acceleration of mass loss from −0.06 ± 0.33 to −1.1 ± 0.06 m yr−1 from 2013, which correlates with the strength of the Pacific Decadal Oscillation. In HMA ice loss is sustained until 2015-2016, with a slight decrease in mass loss from 2016, with some evidence of mass gain locally from 2016-2017 onwards

Acceleration and evolution of faults: An example from the Hunter Mountain-Panamint Valley fault zone, Eastern California

We present new space geodetic data indicating that the present slip rate on the Hunter Mountain–Panamint Valley fault zone in Eastern California (5.0 ± 0.5 mm/yr) is significantly faster than geologic estimates based on fault total offset and inception time. We interpret this discrepancy as evidence for an accelerating fault and propose a new model for fault initiation and evolution. In this model, fault slip rate initially increases with time; hence geologic estimates averaged over the early stages of the fault\u27s activity will tend to underestimate the present-day rate. The model is based on geologic data (total offset and fault initiation time) and geodetic data (present day slip rate). The model assumes a monotonic increase in slip rate with time as the fault matures and straightens. The rate increase follows a simple Rayleigh cumulative distribution. Integrating the rate-time path from fault inception to present-day gives the total fault offset